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I.  PURPOSE  OF  CONTRACT 

The  purpose  of  the  contract  is  to  expose  the  underlying 
phenomena  producing  short  term  frequency  variations  and  to  advance 
the  art  of  short  term  stability  measurement  techniques.  The  re¬ 
sults  of  the  investigation  should  lead  to  the  subsequent  improve¬ 
ment  of  the  short  term  stability  of  high  precision  crystal  oscilla¬ 


tors. 
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II.  ABSTRACT 

Two  definitions  of  frequency  stability  are  discussed.  One 
is  basically  the  variance  of  the  frequency  as  determined  in  an  ob¬ 
servation  time  t.  The  second  definition  involves  the  variance  of 
sample  pairs  of  the  first,  the  sastples  being  separated  by  some  time  T. 
These  definitions  are  shown  to  be  simply  related  to  each  other. 

Since  each  is  derivable  from  the  spectrum  of  the  noise,  it  is  pro¬ 
posed  that  a  convenient  way  of  describing  frequency  stability  is 
to  specify  the  spectrum  of  the  oscillator.  From  the  spectrum  either 
of  the  above,  as  veil  as  other  definitions  of  stability,  may  be 
derived  to  beat  suit  a  particular  application. 

The  remainder  of  the  report  describes  the  work  done  so  far 
on  the  analysis  of  oscillator  non-linearities  as  they  affect  the 
frequency  scatter  introduced  by  noise.  The  spec t rums  of  the  quadra¬ 
ture  components  of  the  noise  for  a  simple  model  of  a  non-linear 
oscillator  have  been  derived.  Further  work  is  required  to  access 
the  importance  of  the  non-linearity  on  the  resultant  frequency 


scatter. 
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III.  PUBLICATIONS,  LECTURES,  REPORTS  AND  CONFERENCES 
Conferences : 

Place:  USASRDL,  Fort  Monmouth,  New  Jersey 
Date:  30  November  1961 

Participants:  Messrs.  Layden,  Schodovsky,  Dr.  Hafner,  (USASRDL) 
Messrs.  Rarity,  Saporta,  Weiss,  (NYU) 

Subject :  Discussion  of  first  quarterly  report  of  subject 

contract. 

Reports : 

Monthly  Progress  Reports 

1.  October  1,  1961  to  October  31,  1961 

2.  December  1,  I96I  to  December  31,  1961 


Lectures : 

Lecture  given  at  New  York  University,  College  of  Engineering, 
Department  of  Electrical  Engineering  Colloquium,  by  Lester 
Saporta  on  December  19,  1961,  entitled,  "Basic  Problems  in 
Frequency  Stability". 
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IV.  FACTUAL  DATA 

A.  Introduction 

In  the  first  quarter  of  this  contract  a  study  of  the 
linear  addition  of  various  noise  spec t rums  to  a  fixed  sinusoid  vas 
made.  The  effects  of  the  noise  on  the  scatter  of  the  observed  fre¬ 
quency  of  the  sinusoid  vas  computed.  The  present  report  extends 
this  work  in  tvo  respects.  First  the  definition  of  frequency  scatter 
has  been  somewhat  generalized  to  permit  the  computation  of  the 
scatter  among  sample  pairs  of  the  observed  frequency  when  these 
samples  are  separated  by  a  given  time.  This  definition  may  be 
more  appropriate  when  a  figure  of  merit  is  desired  in  certain  timing 
applications. 

The  second  extension  of  the  analysis  involves  the  evalua¬ 
tion  of  non-linear  effects  in  the  oscillator  on  the  frequency  scatter. 
This  analysis  is  not  yet  complete  but  the  progress  to  date  is  re¬ 


ported  on. 
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B.  Alternate  Definitions  of  Frequency  Stability 

The  statistical  properties  of  the  phase  accumulated,  in 
a  time  x  by  a  sinusoid  added  to  various  types  of  noise  were  dis¬ 
cussed  in  the  first  quarterly  report  of  this  contract. 

Rie  r.m.s.  value  of  this  accumulated  phase  was  determined 
as  a  function  of  the  measurement  time  t.  This  led  to  a  definition 
of  frequency  stability  as  the  r.m.s.  phase  accumulated  divided  by 
the  time  t. 

In  certain  applications  it  may  be  desirable  to  extend  this 
definition.  For  example  the  velocity  of  an  object  might  be  deter¬ 
mined  by  comparing  the  transmitted  frequency  to  that  of  the  re¬ 
ceived  target  echo.  The  signal  might  be  transmitted  in  bursts  of 
duration  t  corresponding  to  the  time  variable  mentioned  above,  but 
the  time,  T,  between  the  transmission  of  the  burst  and  the  recep¬ 
tion  of  the  echo  might  be  considerably  longer  than  t.  In  this  case 
the  performance  of  the  velocity  measuring  system  depends  on  the 
scatter  among  sample  pairs  of  duration  t  and  spacing  T.  Hie  figure 
of  merit  appropriate  to  the  evaluation  of  an  oscillator  for  this 
purpose  thus  depends  on  both  t  and  T.  This  involves  a  somewhat 
more  complicated  definition  of  frequency  stability  than  the  pre¬ 
viously  used  definition  which  depends  only  on  t. 

The  definition  of  frequency  stability  extended  in  this 
sense  is  easily  derived  from  the  simpler  definition  discussed  in 
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the  previous  report.  This  will  be  developed  in  the  following  para¬ 
graphs. 

In  the  first  quarterly  report  it  was  shown  that  the  effect 
of  adding  noise  to  a  fixed  sinusoid  results,  for  large  signal-to- 
noise  ratios,  is  an  effective  modulation  of  the  phase  of  the  sinu¬ 
soid.  In  the  figure  below,  the  quantity  6( t)  represents  a  typical 
record  of  the  deviation  from  the  nominal  phase  that  would  accrue 


in  the  absence  of  noise.  From  this  diagram, it  may  be  seen  that 


Figure  1  -  Sketch  of  a  representative  record  of  the 
composite  phase  angle  e(t) 

at  time  t,  the  change  in  phase  over  a  short  time  interval  t  is 
Afl(t) .  At  time  t  +  T,  the  change  in  phase  over  the  same  short 
time  interval  t  is  A0(t  +  T).  Ihe  change  in  phase  (or  accumulated 
phase)  may  be  written  in  terms  of  the  absolute  phase  in  the  follow¬ 
ing  way: 

ae(t)  =  e(t  +  t)  -  e(t)  ,  (b-i) 

and:  AS(t  +  T)  =  0(t  +  T  +  t)  -  e(t  +  T)  .  (B-2) 


Corresponding  to  the  change  in  phase  as  given  above,  a  frequency 
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deviation  from  the  nominal  frequency  is  expressible  in  accordance 
with  the  definition  of  frequency  used  in  this  development.  Thus, 


0)(t)  =  -  , 

T 

(B-3) 

and: 

.  .  A6(t  +  T) 

<a(t  +  T)  =  - - - — 

(B-4) 

T 


The  difference  in  frequency  &b(T,  t)  measured  at  the  beginning  and 
at  the  end  of  the  time  interval  T,  may  be  written  as: 

(B-5) 

t)  =  0)(t  +  T)  -  0>(T)  =  -i-[A0(t  +  T)  -  A0(t)] 

Or,  in  terms  of  the  absolute  phase, 

A»(T,  T)  =  -  e3  -  e2  +  e2)  ,  (b-6) 

where,  for  convenience  in  writing  the  absolute  phase,  the  following 
notation  has  been  used: 

ek  -  e(t  +  t  +  t)  , 

e3  =  e(t  +  t)  , 

ez  =  e(t  +  t)  , 


e1  =  e(t) 


(B-7) 
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The  mean  square  value  of  &d(T,  t)  is  therefore: 


a-Ji,  t) 


.  A  2, 


An»(T,  T)  =  —  (ek  -  e3  -  e. 


e1)2 


(B-8) 


Upon  expanding  the  right-hand  side  of  (B-8),  and  noting  that 
0i0j  =  0j0i>  ^'^le  *,°^owinS  expression  is  obtained: 

<T2(T,  t)  = 

0) 


(B-9) 


-±rK 


2  +  «2 


+  ©3  + 


*4  -  20ie2  ‘ 


20163  + 


20^ 


+  202©3  - 


29zek 


-  20304) 


The  quantities  0^  may  be  written  in  terms  of  the  corresponding 
auto -correlation  function  in  the  following  manner: 


02  = 

e|  =  =  el  =  r(o) 

9 

(B-10) 

0ie2 

=  0(t)  0(t  +  t)  =  R(t) 

9 

(B-ll) 

0l©3 

=  0(t)  0(t  +  T)  =  R(T) 

9 

(B-12) 

0le4 

=  0(t)  0(t  +  T  +  t)  =  R(T  +  t) 

9 

(B-13) 

0203 

=  0(t  +  t)  0(t  +  T)  =  R(T  -  t) 

9 

(B-14) 

*2*4 

=  0(t  +  t)  0(t  +  T  +  t)  =  R(T) 

9 

(B-15) 

e3®4 

=  0(t  +  T)  0(t  +  T  +  t)  =  R(x) 

• 

(B-16) 
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Substitution  of  expressions  (B-10)  through  (B-l6)  into  (B-9)  yields: 


(T*(T,  t)  =  — — [4R(0)  -  4R(t)  -  4R(T)  +  2R(T  +  t)  +  2R(T  -  t)]  . 

X2  , 

(B-17) 


Rearranging  expression  (B-17)  results  in: 


<rf(T,  t) 


4(R(0)  -  R(t)]  +  4[R(0)  -  R(T) ] 


-  2[R(0)  -  R(T  +  t ) 3  -  2[R(0)  -  R(T  -  t)] 


(B-18) 


In  the  first  quarterly  report  of  the  contract,  expressions  for  the 
mean  square  value  of  the  accumulated  phase,  for  a  single  sample  of 
duration  t  seconds,  were  derived  for  a  sinusoid  perturbed  by  various 
types  of  random  noise.  This  result  may  be  written  in  terms  of  some 
arbitrary  correlation  time  delay,  5  ,  in  the  following  way: 


<r2(5)  -  2(R(0)  -  R($)]  ,  (B-19) 

AS 

and  may  be  used  to  evaluate  cr^T,  t)  given  in  expression  (B-l8). 

The  procedure  for  carrying  out  this  evaluation  is  to  determine  the 
value  of  <T^J( 5 )  for  f  equal  to  x,  T  and  T  +  r  and  substitute 
these  values  in  (B-l8).  Consequently  expression  (B-l8)  may  be 
written  as: 


t)  -  -i- 


2  <r2(s  ) 

AS 


5  *T 


+  2  <rj(  s ) 

AS 


5  -T 


a-J(S)! 

AS 


-  cr2(S): 

T+T  ^ 


S«T-xJ 


(B-20) 
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2, 

To  illustrate  the  evaluation  of  C^T,  t),  consider  the 

large  signal -to-noise  ratio  approximation  of  the  linear  addition 

2 

of  a  sinusoid  and  random  noise.  The  result  for  when  the 

spectrum  of  the  random  noise  is  shaped  by  a  Gaussian  filter  is 


sketched  in  Figure  2. 


Figure  2  -  Sketch  of  the  variance  of  the  accumu¬ 
lated  phase  for  noise  shaped  by  Gaussian 
filter 


From  the  above,  it  is  seen  how  the  variance  in  the  accumulated  phase 
AS  is  determined  at  the  four  delay  times  t,  T,  and  T  ♦  t.  Know¬ 
ledge  of  the  mean  square  value  R(0)  and  the  auto-correlation  function 

2 

R(S  )  of  the  random  accumulated  phase  is  sufficient  to  determine  <T  (§  )• 


From  the  above  development,  it  has  been  shown  that  an  expression 
for  the  scatter  in  frequency  based  on  sample  pairs  is  derivable 
in  terms  of  the  scatter  based  on  a  single  sample  as  developed  in 
the  first  quarterly  report.  In  view  of  this,  the  scatter  may  be 
determined  from  either  definition  from  a  knowledge  of  the  auto¬ 
correlation  of  the  noise  which  in  turn  is  related  to  the  power 
spectrum  of  the  noise. 
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C.  Non-Linear  Effects  on  Frequency- Scatter 

This  section  of  the  report  is  devoted  to  a  discussion  of 
the  effects  of  noise  In  a  non-linear  oscillator.  The  previous  report 
considered  the  case  of  the  linear  addition  of  various  noise  spectrums 
to  a  fixed  sinusoid.  The  scatter  in  the  accumulated  phase  and  fre¬ 
quency  due  to  the  noise  vas  computed  as  a  function  of  the  observa¬ 
tion  or  measurement  time.  Essentially  the  same  type  of  analysis  is 
required  in  the  non-linear  co.se,  but  the  determination  of  the  effec¬ 
tive  noise  spectrum  is  a  somewhat  more  difficult  problem  than  is  the 
case  for  a  linear  analysis. 

The  following  analysis  represents  the  progress  made  in  this 
problem  to  date.  Further  work  will  be  required  to  permit  the  evalua¬ 
tion  of  the  scatter  in  terms  of  the  frequency  spectrum  which  will  be 
derived. 

The  method  of  approach  was  suggested  by  an  unpublished 
paper  of  Dr.  Eric  Hafner  of  the  United  States  Army  Signed  Research 
and  Development  Laboratories. 

The  circuit  shown  in  Figure  3  represents  the  tank  circuit 
of  a  typicad  oscillator.  The  current  generator  i(t)  provides  for 
the  feedback  due  to  the  presence  of  the  active  element  necessary  to 
sustain  oscillations.  In  general  this  element  vill  be  non-linear 
and  this  is  taken  into  account  by  making  i(t)  a  function  of  the 
instantaneous  voltage  across  the  tank. 
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The  generator  in(t)  represents  the  noise  produced  in  the 
resistor,  R.  For  the  purposes  of  this  analysis  in(t)  will  be  assumed 
to  have  a  white  frequency  spectrum.  For  this  network  we  have 


Ce  +  -2-  ♦  -L 

R  L 


edt  *  i(t)  +  in(t) 


(C-l) 


e  +  GSe  +  co_ 


edt  -  Sf(e)  *  Sin(t) 


(C-2) 


where 


If  the  steady  state  solution  of  the  homogenous  equation  ob¬ 
tained  by  removing  the  noise  source  is  designated  eQ(t),  equation  (C-2) 


becomes 


*0  +  +  ^0 


e0dt  -  Sf(eD)  =  0 


(C-3) 
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Upon  differentiation  vlth  respect  to  time  this  becomes 

*e‘o  +  GSe0  +  e0  -  Sf'(e0)  e0  =  0  ,  (C-4) 


where  the  prime  refers  to  differentiation  with  respect  to  e. 

In  the  weakly  non-linear  case  of  interest  here,  the  steady 
state  solution  will  be  approximately  given  by: 

A 

e0  =  A  cos  u^t  ;  e0  =  -AUq  sin  o^t  ;  *e’Q  =  -Ajj£  cos  u^t 

(C-5) 

Substitution  of  these  results  into  (C-4)  yields 


Ge0  -  f(e0)  eQ  -  0 


(C-6) 


Since  eQ(t)  =  A  cos  u>0t  is  an  even  periodic  function  of  time  with 
fundamental  uJq,  f'(e0)  may  be  expanded  in  a  Fourier  Cosine  Series: 


oo 

V 


where 


f'(e0)  =  Fo  ♦  /  2Fn  cos  nu>ot 
n=l 
p*M> 

Fn  *  -S'  dt 

J-n/(Uo 


(c  -7 ) 


The  differential  equation  (C-6)  must  in  the  steady  state 
be  satisfied  by  each  harmonic  of  the  solution.  Thus  collecting  the 
first  harmonic  terms  after  substituting  (C-5)  and  (C-7)  in  (C-6): 


G*(-AUo  sin  co0t )  -  Fq' (-/hOq  sin  u^t)  +  F2*(-Au0  sin  u^t)  =  0 


or 


G  -  F0  +  F2  =  0 


(C-8) 


15 


Equation  (C-8)  establishes  the  equilibrium  conditions  of 
the  network  without  noise.  When  a  particular  non-linear  function 
f(e)  is  specified,  equation  (C-S)  permits  the  evaluation  of  the 
amplitude  of  the  oscillations. 

Assuming  that  the  steady  state  oscillations  have  been 
established,  the  noise  generator  may  be  reconnected.  Hie  resultant 
solution  for  the  output  voltage  may  be  written  as  the  sum  of  the 
solution  without  noise,  e0(t)  and  some  increment  u(t).  Thus 

e(t)  =  e0(t)  +  u(t)  .  (C-9) 

Substituting  (C-9)  in  equation  (C-2)  results  in 


e<->  +  GSe0  +  as- 


e0dt  -  Sf (eD+u)  +  u  +  GSu  + 


udt  *  Si 


n 


(C-10) 


For  u  very  small 

-  Sf(e0+u)  «  -  Sf(e0)  -  Sf ' (eQ)u  .  (C-ll) 

With  this  approximation,  and  in  view  of  equation  (C-3),  equation  (C-ll) 
reduces  to 


u  +  S[G  -  f'(eQ)]u  +  u>‘ 


udt  =  Si 


n 


(C-12) 


or 


oo 


u  .  Slo  -  ra  2F„  cos  au0t]  u  + 


n«l 


udt  ■  Si 


n 


(C-13) 
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Since  the  tuned  circuit  is  narrow-band,  u  can  have  Important  fre¬ 
quency  components  only  in  the  vicinity  of  u>0.  Neglecting  all  com¬ 
ponents  far  removed  from  aiQj  and  using  the  equilibrium  relation  of 
equation  (C-8),  equation  (C-13)  becomes: 


■fc 

u  -  SFg(l  +  2  cos  2uiQt)  u  +  udt  =  Sin 

V 


(C-14) 


Note  that  equation  (C-l4)  is  a  linear  equation  in  u  with  time  vary¬ 
ing  coefficients.  The  linear  property  permits  the  use  of  super¬ 
position  in  obtaining  a  solution. 

Assume  a  solution  for  u(t)  in  the  form 


Then 


u(t)  =  x  cos  u>0t  -  y  sin  u>0t  (C-15) 

u(t)  =  (x-o^y)  cos  o)0t  -  (y«u>0x)  sin  u>0t  (C-l6) 


The  integral  of  u  with  respect  to  t  is  most  conveniently  obtained 
by  assuming  the  form: 


udt  =  a  cos  u^t  -  b  sin  a^t 


(c-17) 


Differentiating  (C-17)  provides: 


u  =  (a-uiob)  cos  a)0t  -  (iD-Kooa)  sin  u^t 


(C-18) 
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Comparing  (C-l8)  with  (C-15)  results  in 

x  =  a  -  o^b  ;  y  =  fc  +  o^ja  .  (C-19) 

Denoting  by  a  capital  letter  the  Fourier  transform  of  the  corre¬ 
sponding  lo-'er  case  time  function  (C-19)  becomes 

X  =  JcuA  -  <o0B  ;  Y  =  jcuB  +  .  (C-20) 

Equations  (C-20)  may  be  solved  simultaneously  for  A  and  B: 


Juft  +  a>0Y  _  Jo>Y  -  o)0X 

j>  2  *  ?  ? 

<Dq  -  (U  OJq  -  O) 


(c-a) 


Ihe  time  functions  a(t)  and  b(t)  may  be  recovered  by  taking  the  in¬ 
verse  transforms  of  A(u>)  and  B(u>),  i.e.,  a(t)  »  F-1  [A(oj) ] ,  etc. 
However,  since  the  remainder  of  the  analysis  will  be  carried  out 
in  the  frequency  domain  it  is  not  necessary  to  explicitly  perform 
this  step. 

Substituting  (C-15),  (C-l6),  and  (C-17)  into  (C-lk),  col¬ 
lecting  cos  u^t  and  sin  tugt  terms,  and  eliminating  third  harmonic 
components  results  in 


(x  -  <u^y  -  2SFgX  +  cuq  a)  cos  o>0t 


(C-22) 


-  (y  +  u>0x  +  u)q  b)  sin  u>0t  =  Sic  cos  u>Qt  -  Sis  sin  u)Qt  , 


where  the  noise  generator  in(t)  has  been  written  as  the  sum  of 
cos  u^t  and  sin  oj^t  components. 

Separately  equating  cos  aoot  and  sin  u^t  components  on  each 
side  of  the  equation  and  taking  Fourier  transforms,  provides  the 
pair  of  simultaneous  equations: 


O  ■ 2SF*  *  +  1  *  ■  SI= 


(c- 


These  two  equations  may  be  solved  for  X(cu)  and  Y(u>)  in 
terms  of  Ic(o>)  and  I8(<u).  The  resultant  expressions  represent  the 
Fourier  spec t rums  of  the  cosine  and  sine  components  of  the  required 
voltage  perturbation.  The  techniques  for  evaluating  the  mean  square 
values  of  x(t)  and  y(t);  and  the  relation  of  these  quantities  to 
the  amplitude  and  phase  of  the  total  output  voltage  are  similar  to 
those  outlined  in  the  previous  report.  The  actual  integrations 
Involved,  however,  while  straightforward  in  principle  are  unwieldy, 
and  therefore  the  completion  of  this  problem  will  require  some  sim¬ 
plification  of  the  equations  or  perhaps  a  resort  to  graphical  inte¬ 
gration.  The  completion  of  this  phase  of  the  investigation  will 
therefore  be  reserved  for  the  next  quarterly  interval. 
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V.  CONCLUSIONS 

In  discussing  short  term  frequency  stability,  various 
definitions  nay  be  appropriate  to  different  physical  situations, 
lire  relationship  between,  two  definitions  which  appear  to  have  fre¬ 
quent  application  was  explored  in  the  first  section  of  this  report. 
It  was  shown  that  the  statistics  appropriate  to  one  definition  are 
easily  derived  from  those  of  the  other.  This  in  turn  reflects  the 
fact  that  the  statistics  referring  to  either  definition  of  short 
term  stability  are  related  to  the  power  spectrum  of  the  signal. 

In  view  of  this  it  seems  appropriate  to  characterize  an  oscillator 
by  specifying  its  spectrum  rather  than  selecting  any  particular 
definition  of  short  term  stability  as  a  standard.  This  has  the 
advantage  of  permitting  the  user  of  the  information  the  choice  of 
definition  most  appropriate  to  his  own  application.  If  this  ap¬ 
proach  is  used,  the  frequency  scatter  measurements  made  with  the 
multiplier-heterodyne  system  may  oe  regarded  as  a  method  of  deter¬ 
mining  the  auto-correlation  function  of  the  oscillator  frequency 
which  in  turn  may  be  used  to  obtain  the  spectrum.  In  effect  the 
multiplier-heterodyne  system  plays  the  role  of  a  narrow  band 
spectrum  analyzer  using  time  domain  filtering. 

The  second  half  of  this  report  was  devoted  to  a  study 
of  non-linear  effects  in  oscillators  on  the  frequency  scatter 
Introduced  by  noise.  A  model  comparable  to  the  linear  models. 
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which  were  studied  in  the  previous  report,  was  analyzed  following 
a  technique  suggested  by  Dr.  Eric  Hafner,  and  related  to  the 
van  der  Pol  approach.  When  completed,  the  analysis  should  indicate 
whether  a  significant  difference  is  introduced  by  the  non-linearity, 
or  whether  the  simpler  linear  analysis  is  adequate  for  most  pur¬ 
poses.  The  study  has  been  carried  to  a  point  where  the  spectra 
of  the  quadrature  components  of  the  output  sure  available.  Some 
further  work  is  required  to  compare  these  results  with  the  pre¬ 
viously  derived  linear  results. 


VI.  PROGRAM  FOR  THE  NEXT  INTERVAL 


An  important  item  in  the  work  scheduled  for  the  next 
quarter  is  the  completion  of  the  study  of  the  non-linear  effects 
described  in  Section  C.  In  particular  it  is  desirable  to  be  able 
to  indicate  whether  non-linearities  play  a  relatively  important 
or  unimportant  role  in  determining  the  frequency-scatter  due  to 
noise . 

A  second  item  on  the  schedule  calls  for  an  investigation 
of  frequency  determining  techniques  using  a  feedback  system  to 
phase  lock  the  unknown  oscillator  to  a  standard. 

Some  experimental  work  on  the  multiplier -heterodyne 
system  is  also  contemplated.  A  known  frequency  perturbation  will 
be  introduced  into  an  oscillator  and  the  ability  of  the  system  to 
reveal  the  nature  of  the  perturbation  will  be  studied. 


VII.  IDENTIFICATION  OF  KEY  TECHNICAL  PERSONNEL 
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